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The lighter top squark in supersymmetry can live long enough to form hadronic bound 
states if it has no kinematically allowed two-body decays that conserve flavor. In this case, 
scalar stoponium may be observable through its diphoton decay mode at the CERN Large 
Hadron Collider, enabling a uniquely precise measurement of the top-squark mass. The 
viability of the signal depends crucially on the branching ratio to diphotons. We compute 
the ncxt-to-leading order QCD radiative corrections to stoponium annihilation decays to 
hadrons, photons, and Higgs scalar bosons. We find that the effect of these corrections is 
to significantly decrease the predicted branching ratio to the important diphoton channel. 
We also find a greatly improved rcnormalization-scale dependence of the diphoton branching 
ratio prediction. 
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I. INTRODUCTION 

The Minimal Supersymmetric Standard Model (MSSM) [l| with conserved R-parity contains 
a stable lightest supersymmetric particle (LSP). If the LSP is neutral, it could be the cold dark 
matter required by the standard cosmology. The collider signatures of the MSSM generally involve 
missing energy carried away by two LSPs produced in each event. Unfortunately, this suggests that 
there will be no true kinematic mass peaks whose reconstruction would determine superpartner 
masses. In favorable models, it is possible to obtain precision measurements of superpartner mass 
differences and other combinations of masses at hadron colliders by finding kinematic edges from 
decays. However, the overall mass scale of the superpartners will be much harder to obtain precisely 
at the CERN Large Hadron Collider 0,3]. 

In models with a relatively small mass difference between the lighter top squark t± and the 
neutralino LSP Ni, there is an exception that would allow a sharp mass peak. If the lighter 
top squark has no kinematically allowed two-body decays that conserve flavor, then it will form 
hadronic bound states. Among these is stoponium, a stop-anti-stop bound state, which can be 
directly produced at hadron colliders through gluon-gluon fusion. The largest production cross- 
section is for the IS* ( J PC = ++ ) state, denoted in the following by r/j, but other stoponium states 
can contribute to the signal either by prompt decays to the ground state or direct annihilation 
decays. This state will form if 



m t x < m N 1 + m *' 
< + 5 GeV, 



Ci.i) 

[1.2) 



so that the decays t\ — > tN\ and t\ — > bC\ are both kinematically forbidden. (Here C\ stands for 
the lighter chargino mass eigenstate.) These conditions are almost never satisfied in the MSSM 
parameter space with the so-called mSUGRA boundary conditions, but they can easily be satisfied 
in other motivated models. These include "compressed super symmetry" models Q] in which the 
predicted thermal relic density of dark matter is in agreement with that observed by WMAP 
and other experiments due to the enhanced annihilation N\N± — > ti mediated by t-channel 

exchange of top squarks.^ Another class of models consists of those that generate the baryon excess 



12j. In these and other 13J cases, — 



over anti-baryons at the electroweak scale 111 ], 
is necessarily small enough to guarantee the formation of stoponium. 

Stoponium can decay directly by the decays of one of the top-squark constituents through the 3- 
body process t\ — ► WbN±, or if that is kinematically forbidden, through the flavor- violating 2-body 
process t± — ► cN± and/or the 4-body proces s ti — ► ff'bN\. However, the corresponding partial 
widths are many orders of magnitude smaller |l4|-[l8| than the binding energy of stoponium, which 



will be of order a few GeV 19j. Therefore, rjf will decay primarily by annihilation, including the 
possible two-body final states gg, 77, h h°, W + W~, ZZ, Z-y, ti, and bb. Of these, the most 
promising final state, both for detectability over backgrounds and reconstruction of the mass peak, 
is 77, as was first pointed out long ago by Drees and Nojiri 20j, [2l| (See also refs. 22j-[28(] for other 
work related to stoponium at colliders.) The diphoton stoponium signals for both compressed 



The inequalities (jl.ip and (|1.2[1 can also be satisfied in the stop- neutralino co- annihilation region [j|-[lol] of param- 
eter space, but at the price of a much more extreme fine-tuning of input parameters. 
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supersymmetry and supersymmetric electroweak-scale baryogenesis have recently been studied in 
29]. 

In much of the parameter space in which stoponium can form, the gg two-body final state 
leading to hadronic jets dominates. If so, then the leading-order prediction for BR(77) is nearly 
model-independent, and is of order 0.005. The QCD corrections to the bound state annihilation 
decays are quite significant, however, and need to be taken into account along with corrections 
to the production cross-section in order to obtain a realistic estimate of the LHC sensitivity for a 
given model. In this paper, we will calculate the QCD next-to-leading order corrections to S"-wave 
stoponium^ decay into the 77, gg, and h°h° final states. The decay widths into 77 and gg are 
model-independent to leading order, and one may also argue that they are the most important, 
since gg usually dominates the total width, and 77 is the observable signal. However, the h°h° 
final state may dominate in some parts of parameter space, particularly if stoponium is just above 
the threshold for that decay. 



We will proceed following the strategy (and some of the notation) used in ref. 30( where the 
analogous case of quarkonium decay was studied. (The quarkonium annihilation beyond leading 



order was calculated earlier in ref. 3l( , which regulated infrared divergences and mass singularities 
using a gluon mass instead of dimensional regularization.) The S-wave stoponium decay width is 
related to the low-velocity (v — > 0) limit of the stop-anti-stop annihilation cross-section by: 



r( m ^X)=va(i 1 ? 1 ^X)\*(0)\< 



(1.3) 



where ^(0) is the bound-state wavefunction at the origin. [This is often expressed instead in terms 
of the radial wavefunction at the origin, R(0) = y/^r^f(0).] Here, v is the relative velocity of the 
squarks in the center-of-momentum frame. (The same formula (jl.3p holds for excited states with 
angular momentum. Obtaining the decay widths of higher angular momentum stoponium states 
would require keeping contributions at higher order in v.) For diphoton and hadronic final states, 
the cross-section on the right-hand side of eq. (11. 3p is in turn related by the optical theorem to 
the imaginary part of the amplitude for t\t\ — ► t\t\ through two-particle and three-particle cuts. 
For the h°h° final state, we find it easier to just calculate the radiative corrections to the decay 
directly. In both cases, we work in Feynman gauge, and regulate amplitudes using dimensional 
regularization in d = 4 — 2e dimensions. Ultraviolet divergences are indicated separately by writing 
l/euvj while infrared divergences and mass singularities are indicated by 1/eiR, and l/ej R for the 
pole terms. The top-squark propagator is renormalized on-shell, and the QCD gauge coupling will 
be renormalized in the MS scheme. 

An important issue that arises in all calculations of this type is that a(tit* — ► X) obtains 
contributions that are divergent as v — > due to the exchange of massless gluons in diagrams of 
the form shown in figure [TJ The relevant next-to-leading order contribution in QCD is related to 
the leading order contribution by 



A<7«(M* X) 



C F + 0{v") a^ihtt^X) 



(1.4) 



* Although we have in mind the top squark, our results can also be applied to any new strongly interacting funda- 
mental scalar whose width is small enough that it hadronizes before it decays. 
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FIG. 1: Diagrams contributing to tit* — > X with gluon exchange 
between the initial state squarks, leading to Coulomb 1/v singularities 
of the type in eq. (|1.4p . 




where Cp is the quadratic Casimir invariant, 4/3 for SU(3). This Coulomb singularity can be 
absorbed into the definition of the bound state wave- function ^(O). Alternatively, since it is 
universal in character, it cancels when one considers branching ratio observables. This provides a 
useful test of the calculation. 

The rest of this paper is organized as follows. In section [III we find the next-to-leading order 
QCD corrections to stoponium annihilations into gg and 77 final states. Section [TTTl discusses the 
one-loop QCD corrections for the h°h final state. In section ITVl we discuss the numerical impact 
of these results. Section [V] contains some concluding remarks. 



II. DECAYS TO HADRONS AND TO PHOTONS 

In this section, we calculate the next-to-leading order QCD corrections to the gg and 77 partial 
widths. We use the cut method derived from the optical theorem, which allows the direct com- 
putation of the squared amplitude and avoids having to square a matrix element involving many 
terms. To calculate the amplitude for an arbitrary process X — > Y, draw all of the diagrams involv- 
ing X — > X scattering to the desired order, then cut the diagrams through the propagators that 
correspond to the desired final state Y (for tt* — ► gg at leading order, we have the three diagrams 
in figure [3|). These cut propagators are put on mass-shell, using the appropriate Feynman rule 
corresponding to external particles (see Figure [2]). Then the sum of all cut diagrams, multiplied by 
an extra factor of —1 and summed and averaged over spins and colors as appropriate, is denoted 
by M. CVL t and equals the squared amplitude of X — > Y. 

In order to get the partial widths, the contributions to M C ut must be integrated over d- 
dimensional Lorentz-invariant phase space. The individual contribution of a single cut diagram to 
the cross-section is 



1 f.. /Vr 2e d% 

/ V ' ' f 

M cut dLIPSTv- (2.1) 



1 



AEaE b v 



In equation (12. ip . the labels A and B are for the initial-state and / for the final-state momentum 
four-vectors, v is the relative velocity of the initial-state particles, and J dLIPSAr is the integral 
over iV-body Lorentz-invariant phase space. Since we are calculating the annihilation of a bound 
state, we multiply both sides by the relative velocity and set Ea = Eg = mi as the relative 
velocity goes to zero. Therefore, 

vAa = ^5- J M cut dLIPSjy. (2.2) 
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FIG. 2: The Feynman rules used in the calculation of section |TTJ Dashed lines refer to scalar particles, 
solid lines are for fermions, wavy lines are for massless gauge (vector) bosons, and dotted lines are for the 
corresponding ghosts. Dark bars indicate propagator cuts. QED vertices are obtained by the replacement 
ffs-T'fc' 7 — > Qe5 3 k and f abc — > 0. For the four-point scalar vertex, our effective theory neglects Yukawa and 
electroweak couplings that occur in the full MSSM. 



Adding up all of the terms from the appropriate cut diagrams gives the total cross-section multiplied 
by the relative velocity, va, which is related to the partial width in equation (jl.3p . 

At tree-level, the cross section for the annihilation of tit* into a gluon-gluon final state in 
d = 4 — 2e dimensions is 



2m? 

ti 



1^ T(2-e) /V 
T(2-2e) \m 2 



(2.3) 



with ag = (g| / Air) fi~ 2e , where g% is the bare QCD coupling and fi is the regularization mass. The 
diagrams a, 6, and c in figure [3] contribute to this result in the ratio : — 1 : (2 — e). 

The tree-level diphoton cross section can be obtained by the replacement g^T^ 3 — > Qe5 3 k in the 
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FIG. 3: The diagrams whose imaginary parts contribute to the annihilation of t\t\ into gauge bosons at 
leading order. The dark bars indicate where the diagrams have been cut. Diagrams related by permutations 
of external lines and arrow-reversal are not distinguished. 



gg result. For this final state, 

„,.) (¥i ^ ) ^iM« , > (2 . 4) 

m~ 1 (2 — 2e) \ m~ J 

where Q s = +2/3 is the charge of the squark, and a = e 2 /4ir. Therefore, the leading-order decay 
rates are 

r<»> M = f ai ^l, (2.5) 



^rrv, / 128tt 2 I^(0)| 2 

r«»(,^ 77 ) = _^LUL, (2. 6) 



where we have replaced the bare coupling as with the renormalized coupling as, since they are 
equal at leading order. Taking the ratio of these partial widths eliminates the bound state wave- 
function and produces the simple leading-order result 



R( o) = r (0) (m - 77) _ 

~ r(0) (Vi 99) 9a|- ^ () 

The non-vanishing cut diagrams that correspond to the annihilation of t\t\ into gg, ggg, and 
gqq final states at next-to-leading order are given in figure Hf . Many of these diagrams can be cut 
in more than one way. In diagrams with three cut propagators, there is either real gluon emission 
or the pair-production of quarks^, and diagrams with two cut propagators have one-loop integrals. 

In diagrams with three-particle cuts, the principal difficulty is integrating the momentum frac- 
tions of the final-state particles over three-body phase space. To do this, the phase space integrals 
can be reduced to integrals of the form given in the Appendix of ref. [30|. Care must be taken in 



Note that several diagrams not shown in the figure vanish because the color indices of three-gluon final states must 
be antisymmetric by charge conjugation invariance. This is because r\i has C = +1, while a final state with n 
gluons has C = ( — l) n+,lc , where n c is 1 (0) for antisymmetric (symmetric) adjoint color indices [32I ]. 
We include some gqq final state contributions, even though these may be regarded as corrections to qq final states, 
which we do not treat here. The light qq partial widths are suppressed by small Yukawa couplings at leading order, 
and the tt final state is often strongly suppressed by kinematics and couplings. However, the gqq contributions 
from 3-particle cuts in diagrams pi, p2, and p3 in figure [4] cancel large logarithms in the limit of small m q due to 
the gluon vacuum polarization (2-particle cut) contributions from the same diagrams. 




FIG. 4: The diagrams whose imaginary parts contribute to the cross-section for tit* — ► gg at next-to- 
leading order, including ggg and gqq. Each diagram must be cut in all the ways it is possible to put the 
cut propagators simultaneously on-shell, and the propagators that are cut indicate the corresponding two- 
or three-particle final state. Diagrams related by permutations of external states and arrow-reversal are not 
shown, nor are the diagrams with ghost loops that are needed for each diagram with a gauge boson loop. 
We have also not shown several diagrams that vanish for all possible cuts. 
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evaluating diagrams with multiple distinct three-propagator cuts. Diagram d2, for example, has 
two cuts that are not equal. 

Evaluation of the two-particle cuts involves expanding the loop integral from the virtual gluon in 
partial fractions to obtain a set of scalar integrals, which are well-known. A complete set of scalar 
integrals that occur in the calculation can be found in ref. [33] (for a complete set of divergent 



and many finite scalar loop integrals, see ref. 34|). In contrast with the three-particle cuts, the 



phase space integration is quite easy [30|] . Since the cut diagrams do not depend on the final-state 
momentum directions, they are proportional to their contributions to the cross-section 

vAa^ = JL-JWcxtSte), (2.8) 



ti 



where 



1 ( tt \ r(l-e) 

ti , 



*(2) S / dLIPS 2 = - f ^ y . (2.9) 



There is an important simplification that can be made in the massless two-particle cuts of 
diagrams with potential Coulomb singularities (diagrams fl, f2, f3, and f4 in figure Q}. Using the 
identities [30J 

J dLIPS(P = h + k 2 )k^ = ^$(2) (2.10) 
J dLIPS(P = h + k 2 )k1k\ = i L^pfp" _ _1_ P 2^A $(2)> (2 n) 

where k\ and k 2 are the final-state gluon momentum 4- vectors, d = 4— 2e is the number of spacetime 
dimensions, and $(2) the integrated two-body phase space, it is easy to show that, for integrals 
performed in this calculation, dot products of the particle momentum 4-vectors in the numerators 
of loop integrals cannot contain terms linear in the relative velocity v. The effect of the v — > 
divergence comes only from the scalar loop integrals in the calculation, and one may set v = 
everywhere else. 

In Tables HI and [TTl the contribution Cdi ag ram from each diagram to va^(t\i\ — > gg) is given in 
the form 



vAa^ihtl -> gg) = va^ (ht\ -> gg) — /(e) C diagram . (2.12) 

TT 



Here, 




(2.13) 



in keeping with the notation of ref. 30]. Also, Cf is the quadratic Casimir invariant, Ca is the 
Casimir invariant of the adjoint representation, and Tp is the index of the fundamental represen- 
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TABLE I: Results for diagrams not involving propagator corrections to the tree level diagrams. 



tation, given for 577(3) by Cf = 4/3, Ca = 3, and Tp = 1/2. We have combined diagrams ql with 
q2 and q3 with q4 in the table because the individual self-energy diagrams are not proportional to 
the projector g^ v — q^q u /q 2 . 
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TABLE II: Results for diagrams involving propagator corrections to the tree-level diagrams. The function 
h(m 2 j:/m^) is defined in eq. (|2.16p . 



Taking the sum of the diagrams, we find the next-to-leading order result 
(fit? - 99) = - (0) 



m ■> {t x i\ gg) \l + /(e) 
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("light + nt) - 2n t h(m 2 Jm 2 ii ) - - ln(2) j Tp 



(2.14) 



where <5 is either 1 or depending on whether or not the four-point squark interaction in Figure [2] 
is included, § n^ght = 5 is the number of light quarks, and n t = 1 or depending on whether or not 
the top quark is included in the effective theory. In this formula, we have written 



h(r) 



11 4 1 

2 (4 _ r) ^T37 _ 8 _ | i n (i + + 2 ln(2) 



(2.15) 
(2.16) 
(2.17) 



2 / 2 

r = m t /m^. 

The function /i(r) is defined so that it parametrizes the effects of a non-zero top-quark mass. In 



In the MSSM, 5 = 1. However, one can imagine non-supersymmetric theories with fundamental strongly interacting 
scalars, in which these formulas would apply with 5 — 0. 
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the limit that the top quark is massless compared to the top squark we have /i(0) = 0, and when 
the masses are identical we have h(l) = — | + | ln(2). 

The one-loop order correction to the diphoton cross-section can now be found simply by dropping 
the diagrams that involve gluon self-coupling or real gluon emission (equivalent to setting Ca = 0) 
as well as any vacuum polarization diagrams (which no longer involve strong couplings), then 
making the replacement g^Ti — > QeS J k at vertices to change gluons into photons. Following this 
procedure, we find 



va 



(i) 



(Mi -» 77) = va^ (htl -» 77) 



l + f(e)^C F 

IT 



---\5 
2 



7T 2 11 7T 

+ 2 In 2) 

v 2 8 y J 



(2.18) 



In the MS renormalization scheme, the bare coupling as is written in terms of the renormalized 
running coupling as(Q) using 



as = as 



47r \euv 



+ In (4V/Q 2 ) " IE 



(2.19) 



The gluon cross section as a function of the renormalized MS coupling as and the renormalization 
scale Q can therefore be written as 



va W (ifi - gg) = va^> (fctj - gg) { 1 + ^ 



a s 



2 V4m? 



.'vr 2 7 ^ 2 (I vr 2 \ . , 



199 IStt 2 . , 
' 18" ~ ~W ' A 



+ (-y(niight+^)-2n 4 / l (m 2 /m|)-iln(2)^T F j. (2.20) 



For the diphoton final state, one can simply replace the bare coupling 2g by the MS coupling as 
in eq. (|2.18p to obtain the corresponding renormalized result for va^ (t\tl — > 77), since the QCD 
coupling does not appear in the tree-level result in that case. 

The ratio of partial widths is now obtained at next-to-leading order by 



r (1) fa^77) 

r^) (rft — > hadrons) 



(2.21) 



(Here we write "hadrons" to subsume the gg, ggg, and the partial gqq parton-level contributions.) 
The 1/v Coulomb singularity does not appear in this ratio, since it can be absorbed into a redef- 



inition of the bound-state wavefunction factor 37J, at least at order as in the approximation of 
a Coulombic bound state, and the redefined bound-state factor in turn cancels from the ratio of 
decay rates. This redefinition simply removes the 1/v part as it is expanded to next-to- leading 
order in as- There may remain some small residual dependence on 1/v proportional to a s , due to 
the fact that the bound-state potential is actually not exactly Coulombic, but this is beyond the 
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FIG. 5: The tree-level diagrams for the annihilation oit\t\ into h°h°. 
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(2.22) 



An interesting feature of this result is that the term proportional to Ca is the same for the stopo- 
nium and the corresponding quarkonium calculation as can be seen by comparing^ eq. (I2.22|) to 
eqs. (4.1)-(4.4) of ref. 0. 



III. DECAYS TO HIGGS SCALAR BOSONS 

We now consider the one- loop radiative corrections to m — * h°h°, where h° is the lightest Higgs 
scalar boson in supersymmetry. In this case we calculate the partial widths directly rather than 
using the cut method. The tree-level diagrams contributing to this annihilation decay are shown 
in figure EQ T he corresponding annihilation cross-section in the v — * limit can be written as 21] 
(see also HQ]): 



N r . 



647rm? 



(l-m 2 h0 /ml) 1 / 2 {\ 1 + \ 2 + \ 3 ) 2 , 



where, in the notation of ref. [3g |. the effective couplings are: 
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(3.4) 



' Note that in ref. [3C| is the reciprocal of our definition. 
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FIG. 6: The one-loop QCD diagrams for the annihilation of tit* into h°h° 
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^ + £ - 1 - ln(l - h/2) - 2=nh(h) - \k 2 {h) 


H3d 


_J_ + 1 + i n(2 ) _ ^ ln(2 - h ) + ±h(h) + k 2 (h) 



TABLE III: Results for one-loop radiative corrections to — ► h°h°, corresponding to the diagrams in 
figure [6] and appearing in eqs. (|3.7|) , (|3.8|) . 



Using eq. I2.6p . the tree-level ratio of 77 to partial widths is therefore 

2048vr 2 a 2 



r(°)(77)/r(°)(/ l / l ) 



81(Ai +X2 + X3)' 



1 



m h0 /m k ) 



2 N-l/2 



(3.5) 



The one-loop QCD corrections to the stoponium decay to are due to the diagrams shown 
in figure EJ Note that there are no gluon emission diagrams to this process at this order, because 
the initial and final states are both color singlets. In this paper, we will neglect^ in A3 the con- 
tribution of the heavier top-squark mass eigenstate t 2 , which would otherwise entail a somewhat 
more complicated kinematic loop integration. At one-loop order in QCD and in the limit of small 
v, the corresponding t±t* — ► h°h° result is: 



h°h°) 



N r . 



QAirm 2 

ti 



(1 - ml»/m 2 h ) l/2 [c(Ai + A 2 ) + C% 



(3.6) 



The decay of stoponium to h°h° is especially important in electroweak-scale baryogenesis models [ill], 0, [3 that 
require light t%. In those models, is necessarily very large, making this approximation extremely good. 
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where, in terms of the individual diagram contributions given in Table Hill 



a s 



C = 1 + C F — /(e) (Cm« + Cjnt) , 

7T 

C" = 1 + Cf— /(e) (C#3a + C#3b + Ch3c + Cnsd) 

IT 



(3.7) 
(3-8) 



The couplings Ai,2,3 are given by the same formulas as eq. (|3.2p - (|3.4|) . but with bare couplings, 
denoted by A^o^f*, A^o^*, and X^o h o h o, in place of their unhatted counterparts. Also, 



h = mlo/m^, 



(3.9) 



and we have defined functions 



k x {h) 
k 2 (h) 



Bt&n~ 1 {hB/(2 - h)), 
2-h 
~2A 



Re 



Li, 



l + A 



I- A 



+-Li 2 (1 



1 + iABJ \l + iAB 

2A \ 1 / 2A 
-Li 2 1 + 



2-h 



2-h 



(3.10) 



(3.11) 



with A = \/l — h and B = y/A/h — 1, and Li 2 (x) is the dilogarithm function (also known as the 
Spence function). These functions have values k\(0) = 2 and and £; 2 (0) = — 7r 2 /8 for the extreme 
limit m h o <C m^, and k\(l) = n/V3 and fc 2 (l) = — 7r/2\/3 near threshold for the decay. The 
calculation of these diagrams relies on loop integrals that can be found from refs. 33], |36J] by taking 
the v —* limit with appropriate special cases of momenta and masses. Note that the 1/eiR, poles 
cancel in eqs. (|3.7p and (|3.8p . as required. 

The bare couplings can be written in terms of the renormalized running MS scheme couplings 
^hPhPiiil' ^<f>°tit*> anc ^ / V'h°/i ' a ^ one-loop order in QCD, as: 



A 



A </>°tit* 

A^OfeO/jO 



A 



hOhOtJl 



A^OfeOftO, 



_3«s/J_ + ln( 4 2 /Q2) _ 7£ 
47r Veuv 

+ln(4V/Q 2 )-7^ ^ 

47T Veuv 



(3.12) 

(3.13) 
(3.14) 



which eliminates the 1/exjv dependence of the result up to terms of 0(a 2 s ). It follows that 



va^ihtt -> W) 



iV, 



647rm? 



(1 - m^/m?) 1 / 2 ^ (Ax + A 2 + K'A 3 ) 2 , 



(3.15) 



where Ai j2j 3 are as given in eqs. (|3.2p - (|3.4p . with A^o^f*, A^o^f*, and A^o/jO/jO taken to be MS 
couplings, the masses renormalized on-shell, and 



K 



l + C F 



as 

7T 



7T 



+ - ln(Q 2 /4m 2 ) - 1 + 3m(2) 



(3.16) 
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K' = 1 + C F 



as 

7T 



ln(g 2 /4m|) + | - \ ln(2) - 21n(l - ty2) 



1 



(3.17) 



By comparing eqs. (|2.18p . (|3. 1 1) . and (|3.15p . we obtain 

rW(77) _ rW(7 7 ) (A 1 + A 2 + A 3 ) 2 / a* 
rW (W) r(°) ( W) (Ai + A 2 + £T'A 3 ) 2 V F vr 



ln(Q 2 /4m?) + 4 + 51n(2) ), (3.18) 



written in terms of MS couplings and on-shell masses. 



IV. NUMERICAL RESULTS 



In this section, we will examine the numerical impact of the radiative corrections found above in 
typical realistic cases. In many models, the dominant final state is hadrons coming from the gg or 
ggg parton-level process at next-to-leading order in QCD. As found in ref. 29j], this often amounts 
to 90% or more of the total decay width. (In particular, even when it is kinematically allowed, the 
it final state is typically dominated by the gg final state and has a branching ratio of only a few 
per cent, unless there is a resonant Higgs exchange contribution.) Therefore, for simplicity we will 
begin by considering the idealized case that only diphoton and gluon-induced hadronic final states 
are included. 

The relevant QCD group theory invariants for SU(3) C are Ca = N c = 3, Cf = (A" 2 — X)/2N C = 
4/3, Tp = 1/2, and nu g ht = 5. At two-loop order, the QCD coupling runs according to [38|] 



Q 



das 
~d~Q 



bo 2 
"2^ 



23 2 

T ~ 3 
58 19 



h 3 

4^' 
1 



nt 



6 



11 



rif n, 

3 3 6 



(4.1) 
(4.2) 
(4.3) 



where nt and n s are either or 1 depending on whether or not the top and stop, respectively, are 
included in the effective theory. Although we have not included the QED correction diagrams to 
the annihilation process, we do incorporate the running of the QED coupling according to 



a{Q) 



a(Qo) 



uEM 



1 + (b§ M a(Qo)/2vr) HQ/Qo) ' 



80 16 4 

"T; 7T n t - 7i n s 



(4.4) 



We take as inputs M z = 91.18 GeV, m t = 172.5 GeV, and Og ] {M z ) = 0.118, a^(M z ) = 1/128.0 
as the running MS coupling inputs in the 5-quark Standard Model. We then determine MS cou- 
plings as(iTif ) and a(m^ ) by running below Q = mt using the 5-quark two-loop renormalization 
group equations, and, if > mt, running above Q = mt using the 6-quark renormalization group 
equations. 
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FIG. 7: The dependence of the ratio of decay widths of stoponium into 77 + X and into hadrons is shown 
at tree level (R^°\ dashed lines) and at next-to-leading order in QCD, (i?^ 1 - 1 , solid lines), as a function of 
the MS renormalization scale Q. The left panel shows the results for mi = 120 GeV and the right panel 

238 GeV and 447 GeV, respectively. 



for mi i = 225 GeV, corresponding to approximately m rji 



Having determined the input parameter values at Q = , we then match ^ onto and work in 
the effective theory in which the lighter top squark is always included (even for Q < m^ ), so that 
as{Q) and a(Q) are obtained by running the renormalization group equations eq. (|4.ip and (|4.4j) 
with n s = 1 and either i%t = 1 (if > mt) or m = (if m t - < mt). We then evaluate the ratios 
of decay widths of stoponium into 77 + X and into hadrons at tree-level and at next-to-leading 
order, obtained from eqs. (|2.7p and (|2.22j) . respectively, self-consistently using the same value for 
n t . 

In figure we compare the predicted ratios and as a function of the choice of MS 
renormalization scale Q, for two different choices of stop mass 120 and 225 GeV, corresponding 
to stoponium masses m rji = 238 and 447 GeV. The leading-order prediction has a strong scale 
dependence, which is not surprising since it is inversely proportional to o? s . The next-to- leading 
order prediction for the branching ratio to photons is considerably smaller for all values of Q 
considered. For the choice Q = mi, the decrease is roughly 30%. This may be considered 
somewhat unfortunate for the observability of the signal. We also note that the Q dependence of 
the predicted branching ratio is much improved by the next-to-leading order calculation. 

Figure shows the ratio of diphoton and hadronic decay widths as a function of the stoponium 
mass, computed at leading order and at next-to-leading order. The results are shown for two choices 
of the renormalization scale Q = /2 and Q = 2m^ . Again we see that the diphoton branching 
ratio predicted by the next-to-leading order calculation is smaller than predicted at leading order. * 
Taking into account the fact that the contributions of other final states (h h , W + W~ , ZZ, and 
ti) are quite model-dependent but can only decrease the diphoton branching ratio, we conclude 



We neglect threshold corrections when matching between different effective field theories, because they are of 
relative order a| and numerically small compared to other sources of error and uncertainty. 

Small kinks are visible in the prediction curves in Figure [8] at m r = 2mt, due to our use of the 5-quark (6-quark) 
effective theory below (above) that threshold. 
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that BR( 7 7) is not more than about 0.25% for m Vi = 200 GeV and 0.35% for m n{ = 500 GeV. 

It is somewhat more difficult to make a model-independent statement about the impact of 
the QCD radiative corrections for the h h final state. However, one can note that the factor 
K' appearing in eq. (|3.15p is generally quite close to unity. Typical numerical values found from 



eq. ([ST7D are 0.006 < K' - 1< 0.021 for Q 



m- tl and 0.050 < K' 



1 < 0.065 for Q = 2m~ t . The 



lower range in each case occurs closer to threshold for the decay. Neglecting this small deviation 
from unity, one finds from eq. (|3,18p 



BR( 77 ) 



BR(h°h°) 



BR( 77 ) 



NLO 



BR(h°h°) 



LO 



as r3 



C F — 

7T 



-ln(Q74mf i )+4 + 51n(2) 



(4.5) 



yielding a correction of order (—23%, —32%) for Q/m^ = (1,2), respectively. This correction 
is comparable to, and has the same sign as, the ratio of the diphoton to the hadronic branching 
ratios. However, we note that the BR(/i°/i°) is itself highly mo del- dependent. It is generally quite 
small in the compressed supersymmetry models where top squarks mediate the annihilation of dark 
matter in the early universe, but can be very important in the models with light stops motivated 
by baryogenesis. 



V. OUTLOOK 

In this paper, we have calculated the stoponium decay rates to gg, 77 and h°h° at next-to- 
leading order in QCD, in order to obtain the branching ratio to two photons. Our calculation 
applies to the S-wave J PC = ++ states, including the ground state that has the largest direct 
production cross-section. We have not included the W + W~ and Z°Z° partial widths, which are 
often the next-largest contributions to the total stoponium decay rate, and have neglected the 
remaining decay channels, which have branching ratios of no more than a few percent at leading 
order. t However, it is worth noting that we can approximate the radiative corrections to the 



This includes the tt final state, which is dominated by the gg final state and so has a branching ratio of only a 
few percent even when is much larger than mt, unless there is a resonance Higgs exchange contribution. 
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decay widths to W + W~ and Z°Z°, provided that the contributions from t-channel b and t squark 
exchange can be neglected. In this approximation, the radiative corrections come exclusively from 
4-point and s-channel diagrams similar to Hla through H2b of figured! and it follows that the ratio 
of partial widths takes exactly the same form as eq. (j4.5h with the Higgs branching ratio BR(/i°/i°) 
replaced by BR(W /r+ W^~) or BR(Z° Z°) as appropriate. 

This simplification applies to both of the model scenarios that motivate our study. In order 
to generate a relic density of dark matter in agreement with observation, models with compressed 
supersymmetry should have an LSP between about 170 and 270 GeV and a lightest stop mass 
eigenstate approximately 25 to 100 GeV heavier, which implies a stoponium mass very roughly 
between 400 and 700 GeV 0]. Within this range, there is extra suppression in the W + W~ and 
Z°Z° decay channels due to a cancellation between the terms of eq. (A. 9) in ref. 29[] that come 
from the 4-point and s-channel Higgs exchange diagrams (see also figure 5 of that paper for typical 
branching ratios as a function of stoponium mass in compressed supersymmetry). Even with a 
stoponium mass as high as 600 GeV, the combined and Z°Z° decay rates represent no 

more than ten percent of the total width. In any case, the radiative corrections to the model- 
dependent terms arising from t-channel squark exchange are never large enough to make a big 
difference in the total BR(77). 

In contrast, in models that feature electroweak-scale baryogenesis, the combined contribution 
to the stoponium d ecay rate from W + W~ and Z°Z° can be as large as about 50% at leading order 
(see figure 8 in ref. 29(]). However, in this scenario W + W~ and decays from t-channel squark 

exchange are highly suppressed because the light stop consists almost entirely of the right-handed 
gauge eigenstate and the mass of the lightest sbottom is much greater than the mass of the lightest 
stop (for leading order amplitudes, see eqs. (A. 5) through (A. 10) of ref. [29[]). Therefore the ratios 
of the 77 partial width to the W + W~ and Z°Z° partial widths will take the same form as eq. 
P~5]h with BR(W + W~)ot BR(Z°Z°) in place of BR(/i°/i°), to a very good approximation. 

It should be noted 2l| that the S'-wave excited states of stoponium can be produced at the 
LHC with either direct decays to two photons or decays (with emission of soft mesons or photons) 
to lower stoponium states before annihilation. The treatment of direct annihilation decays of the 
excited S-wave states is the same as for the ground state. These contributions will be essentially 
indistinguishable in the 77 invariant mass signal, since the binding energy differences between such 
states will be less than the experimental resolution. However, part of the resulting signal will be 
lost when the excited S'-wave states decay to P-wave states that decay directly [21(. Evaluating 
the contributions from excited states will require a more detailed understanding of the stoponium 
spectroscopy. 

Although radiative corrections decrease the branching ratio significantly, stoponium annihila- 
tion to 77 is still a viable signal at the LHC. The approximately 30-35% reduction in the diphoton 
branching ratio found in this paper is likely to be offset by an enhancement due to radiative correc- 
tions to the cross-section for stoponium production in hadron colliders; note that the corresponding 
calculation for toponium production results in a K factor of roughly 1.3 or more 39f] . We plan to 
report on the corresponding results for stoponium in a future paper. In any case, with sufficient 
integrated luminosity, stoponium should be visible, and the diphoton decay channel remains a 
unique opportunity for the direct and precise measurement of superpartner masses. Moreover, a 
measurement of the rate for pp — > — > 77, taking into account the radiative corrections to both 
production and decay, will provide interesting and useful information about the stoponium system. 
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